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ABSTRACT 

We present orbital parameters for six double-lined spectroscopic binaries 
(i Pegasi, oj Draconis, 12 Bootis, V1143 Cygni, j3 Aurigae, and Mizar A) and 
two double-lined triple star systems (k Pegasi and rj Virginis). The orbital fits 
are based upon high-precision radial velocity observations made with a dispersed 
Fourier Transform Spectrograph, or dFTS, a new instrument which combines in- 
terferometric and dispersive elements. For some of the double-lined binaries with 
known inclination angles, the quality of our RV data permits us to determine 
the masses Mi and M 2 of the stellar components with relative errors as small as 
0.2%. 

Subject headings: binaries: spectroscopic, techniques: radial velocities, instrumenta- 
tion: spectrographs 



-3 - 



Introduction 



For the past several years, our research group has been developing a new optical 
spectrograph concept called the dispersed Fourier Transform Spectrograph, or dFTS. 
The instrument design merges a traditional Fourier Transform Spectrometer (FTS) with 
a dispersive grating spectrograph, such that the interferometer output is divided into 
thousands of narrowband channels, all operating in parallel. This multiplex advantage 
boosts the effective throughput of the system by a large factor, making the dFTS 
competitive with echelle spectrographs for spectroscopic analysis of stars, particularly 
measurement of their radial velocities (RVs). 



Hajian et al. 



(120071 ) describes our prototype device, dFTSl, and explains the underlying 
theory and hardware implementation in detail. Based upon our commissioning observations 
with dFTSl, we subsequently designed and built a second-generation version, dFTS2, which 
we deployed t o the Stewa r d Ob servatory 2.3-meter Bok Telescope for a year-long observing 



campaign. In 



BehretaL 



(120091 ). we discuss the dFTS2 hardware and present velocimetry 



measurements of RV standard stars and single-lined spectroscopic binary stars (SBls). 

In this paper, we describe the results from our dFTS2 observations of double-lined 
spectroscopic binaries (SB2s) and double-lined triple systems. SB2s provide one of the 
best means for measuring the masses of stars: given an accurate RV curve for each stellar 
component and the inclination angle i of the orbital plane to the observer's line of sight, 
we can derive the component masses using Kepler's Third Law. Traditional spectroscopic 
observations with an echelle spectrograph and thorium-argon calibratio n source can achiev e 

)w-lined stars ( 



Skulian et al. 


2004 


precision, 


Konacki 



Tomkin fc Fekel 



(12005 



2006 



Fekel et al. 



2007 



Ramm 



Ramm et al. 



2004; 



20081 ) . For greater 



20091 ) has developed a technique using an iodine absorption cell, 



with which the RVs of a spectroscopic binary can be measured at the 0.005-0.010 km/s 
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level. Our dFTS2 instrument, in contrast, achieves high RV p r ecisio n and stability without 



BehretaL 



(120091 ) . we measure the RVs 



a superposed reference spectrum. As described in 
of non-binary stars and single-lined binaries to 0.01-0.03 km/s, and anticipate even better 
performance once thermal stability issues in our instrument design have been addressed. 



2. Data acquisition and RV analysis procedure 

The data reported in this paper were collected between October 2007 and June 
2008 during bright-time observing runs at the 2.3-meter Bok Telescope of the Steward 
Observatory on Kitt Peak. An observation of a given SB2 target consisted of 500 
exposures spanning a range of interferometer delays corresponding to a spectral resolution 
of approximately 50,000. Each exposure lasted 1.0 to 4.0 seconds in duration, depending 
on the star's brightness and the atmospheric seeing and opacity. Each scan also required 
a total overhead time of approximately four minutes, independent of exposure time, for 
CCD readout and moving to the next delay position. Our targets, listed in Table [TJ were 
chosen because of their relative brightness and short periods, so that we could acquire many 
observations per star during the limited duration of this initial observing campaign. 

To measure double-lined RVs from our interferogram dat a, we employed a varia tion 



1994T ). 



of the standard two-dimensional cross- correlation technique ( IMazeh &: Zucker 
Instead of transforming our interferograms into spectra, and then cross- correlating 
model templates against each observed spectrum, we convert the template spectra into 
template interferograms, and then compare those model interferograms to the observed 
interferograms. Because interferograms add linearly, we can compute sequences of 
single-lined interferograms for the A and B components of a binary, spanning a range 
of radial velocities for each, and then add together A and B interferograms to create a 
two-dimensional grid of double-lined model interferograms. Calculating the \ 2 difference 
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between the models and the observed interferogram data, we construct a map of fit quality, 
where the minimum point indicates the best-fit solution for component velocities V\ and v 2 , 
and the projection of the x 2 — Xmm + 2-30 contour line onto the V\ and v 2 axes provides 1-cr 
error bars on each radial velocity point. The shapes of the \ 2 contours were elliptical with 
minor and major axes aligned to the v% and v 2 axes, indicating no significant covariance 
between the component velocity errors. 

For optimal results, the template spectra must be well matched to the ac- 
tual spectra of the two stellar components. We generated synthetic sp ectra us- 



ing the SPECTRUM spectral synthesis package ( IGray &: Corbally 



19M) (see also 



http://www.phys.appstate.edu/spectrum/spectrum.html), and then varied the transi 



tion strength log gf of each line and the projected rotation velocity v sin i for each stellar 
component to minimize the x 2 difference between the model and the observed data. A 
final template spectrum for each component was calculated from a median-filtered average 
of transition strengths from the individual observations, and the final two-dimensional 
RV cross-correlation was then performed using these templates. We found that we could 
derive relatively precise and self-consistent v sin i values from our interferograms (as listed 
in Tabled]); a future paper will explore the use of dFTS data to measure stellar rotation 
velocities and other line broadening mechanisms. 



-6 - 



Table 1. Spectroscopic binary targets observed with dFTS2. 



star 


V magnitude 


spectral type 


P (days) 


v sin 1^4 g (km/s) 


references 


i Peg 


3.8 


F5V 


10.21 


7.6, 7.2 


Fekel & Tomkin (1983) 


lj Dra 


4.8 


F5V 


5.28 


7.1, 6.6 


Mavor & Mazeh (1987), Fekel et al. (2009) 


12 Boo 


4.8 


F8IV 


9.60 


13.1, 10.4 


Boden et al. (2005), Tomkin & Fekel (2006) 


VI 143 Cyg 


5.9 


F5V 


7.64 


23.6, 36.2 


Andersen et al. (1987) 


Aur 


1.9 


A2IV 


3.96 


34.5, 35.0 


Smith (1948), Pourbaix (2000) 


Mizar A 


2.3 


A2V 


20.54 


32.6, 36.2 


Fehrenbach & Prevot (1961). Pourbaix (2000) 


k Peg 


1.2 


F5IV 


5.97 


7.1, 47.9 


Mavor & Mazeh (1987), Haiian et al. (2007) 


rj Vir 


3.9 


A2V 


71.79 


5.1, 4.5 


Hartkopf et al. (1992). Hummel et al. (2003) 
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Table 2. Radial velocity data measured with dFTS2. 



star 


HJD - 2,400,000 


phase 


Vi (km s- 1 ) 


Vi error (km s 1 ) 


V 2 (km s" 1 ) 


V2 error (km s - 1 ) 


i Peg 


54,400.6388 


0.109 


33.302 


0.017 


-64.348 


0.080 


i Peg 


54,401.7122 


0.214 


6.584 


0.019 


-21.902 


0.091 


t Peg 


54,401.7269 


0.215 


6.235 


0.017 


-20.958 


0.084 


i Peg 


54,402.6529 


0.306 


-20.809 


0.021 


22.781 


0.109 


t Peg 


54,403.6604 


0.404 


-44.161 


0.020 


59.933 


0.103 


i Peg 


54,404.6147 


0.498 


-52.709 


0.032 


73.279 


0.196 


i Peg 


54,604.9660 


0.115 


31.993 


0.017 


-62.573 


0.080 


i Peg 


54,606.9635 


0.311 


-22.230 


0.015 


24.831 


0.067 


t Peg 


54,634.8758 


0.044 


42.325 


0.015 


-78.934 


0.067 


i Peg 


54,634.8959 


0.046 


42.119 


0.016 


-78.687 


0.072 


i Peg 


54,635.9132 


0.145 


25.341 


0.017 


-51.490 


0.076 


t Peg 


54,635.9342 


0.147 


24.823 


0.016 


-50.598 


0.075 


i Peg 


54,637.8866 


0.339 


-29.723 


0.017 


36.656 


0.073 


i Peg 


54,637.9068 


0.341 


-30.242 


0.017 


37.434 


0.080 


i Peg 


54,638.9426 


0.442 


-49.514 


0.015 


68.427 


0.066 


i Peg 


54,638.9628 


0.444 


-49.713 


0.015 


68.752 


0.064 


oj Dra 


54,547.9814 


0.160 


-48.332 


0.027 


29.389 


0.049 


oj Dra 


54,549.9958 


0.542 


19.025 


0.026 


-53.678 


0.046 


u> Dra 


54,578.9388 


0.024 


-44.148 


0.031 


24.305 


0.056 


oj Dra 


54,579.9848 


0.222 


-41.842 


0.027 


21.495 


0.050 


lu Dra 


54,603.8020 


0.733 


13.203 


0.024 


-46.284 


0.045 


6ti Dra 


54,603.9267 


0.756 


9.247 


0.029 


-41.483 


0.054 


oj Dra 


54,604.9389 


0.948 


-31.906 


0.026 


9.346 


0.049 


ui Dra 


54,605.8126 


0.114 


-49.790 


0.023 


31.110 


0.041 


a; Dra 


54,605.9378 


0.137 


-49.461 


0.026 


30.731 


0.048 


u> Dra 


54,634.9202 


0.627 


22.606 


0.025 


-57.902 


0.046 


a) Dra 


54,634.9758 


0.637 


22.298 


0.023 


-57.528 


0.040 


oj Dra 


54,635.6924 


0.773 


6.250 


0.022 


-37.884 


0.039 


ui Dra 


54,636.6824 


0.960 


-34.305 


0.024 


12.288 


0.042 


oj Dra 


54,636.9410 


0.009 


-42.260 


0.022 


22.042 


0.039 


u> Dra 


54,637.6927 


0.152 


-48.736 


0.027 


29.976 


0.049 


oj Dra 


54,637.9312 


0.197 


-45.017 


0.026 


25.330 


0.047 
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Table 2 — Continued 



star 


HJD - 2,400,000 


phase 


Vi (km s- 1 ) 


Vl error (km s x ) 


V 2 (km s- 1 ) 


V2 error (km s x ) 




12 Boo 


54,547.9153 


0.591 


-29.048 


0.065 


50.273 


0.072 


12 Boo 


54,548.9127 


0.694 


-48.375 


0.076 


70.324 


0.083 


12 Boo 


54,549.9111 


0.798 


-51.801 


0.080 


73.760 


0.093 


12 Boo 


54,577.7996 


0.702 


-49.395 


0.049 


71.239 


0.056 


12 Boo 


54,579.8185 


0.912 


-18.824 


0.060 


39.831 


0.070 


12 Boo 


54,579.8722 


0.918 


-16.030 


0.054 


37.061 


0.062 


12 Boo 


54,602.8623 


0.312 


50.008 


0.070 


-31.349 


0.079 


12 Boo 


54,603.7304 


0.402 


22.867 


0.073 


-3.201 


0.084 


12 Boo 


54,604.7012 


0.503 


-6.315 


0.061 


26.942 


0.071 


12 Boo 


54,605.7012 


0.607 


-32.713 


0.055 


53.903 


0.065 


12 Boo 


54,606.6579 


0.707 


-49.876 


0.063 


71.703 


0.071 


12 Boo 


54,635.7226 


0.733 


-52.290 


0.064 


74.099 


0.073 


12 Boo 


54,638.7263 


0.046 


58.158 


0.084 


-39.386 


0.092 


V1143 Cyg 


54,603.8512 


0.742 


6.959 


0.260 


-39.988 


0.490 


V1143 Cyg 


54,604.8585 


0.873 


57.691 


0.200 


-92.635 


0.333 


V1143 Cyg 


54,605.8830 


0.007 


59.725 


0.191 


-95.517 


0.309 


V1143 Cyg 


54,606.8736 


0.137 


-69.958 


0.250 


37.136 


0.445 


V1143 Cyg 


54,634.8268 


0.796 


22.287 


0.168 


-55.975 


0.320 


V1143 Cyg 


54,636.8467 


0.060 


-29.744 


0.205 


-0.671 


0.339 


V1143 Cyg 


54,638.8573 


0.323 


-65.668 


0.192 


33.420 


0.335 


Aur 


54,400.9449 


0.134 


54.375 


0.151 


-91.605 


0.149 


Aur 


54,401.9169 


0.379 


-96.561 


0.158 


63.013 


0.156 


Aur 


54,402.9819 


0.648 


-81.719 


0.193 


48.833 


0.195 


Aur 


54,402.9909 


0.651 


-81.088 


0.169 


47.313 


0.170 


Aur 


54,403.8922 


0.878 


60.619 


0.289 


-98.087 


0.286 


Aur 


54,403.8982 


0.880 


61.415 


0.203 


-98.274 


0.196 


Aur 


54,404.8387 


0.117 


62.259 


0.292 


-99.681 


0.290 


Aur 


54,404.8542 


0.121 


60.690 


0.241 


-97.396 


0.236 


Aur 


54,404.8680 


0.125 


58.476 


0.253 


-95.754 


0.244 


Aur 


54,487.7015 


0.042 


86.798 


0.292 


-124.300 


0.274 


Aur 


54,487.7106 


0.044 


86.270 


0.237 


-124.136 


0.236 


Aur 


54,487.7191 


0.046 


86.305 


0.381 


-123.834 


0.374 
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Table 2 — Continued 



star 


HJD - 2,400,000 


phase 


Vi (km s- 1 ) 


Vi error (km s 1 ) 


V 2 (km s- 1 ) 


Vi error (km s v ) 


Aur 


54,487.7530 


0.055 


84.279 


0.550 


-122.404 


0.556 


Aur 


54,488.7643 


0.310 


-56.848 


0.217 


23.889 


0.209 


Aur 


54,488.7731 


0.312 


-58.307 


0.206 


24.076 


0.203 


Aur 


54,488.7818 


0.315 


-59.895 


0.200 


26.191 


0.196 


Aur 


54,491.6686 


0.044 


85.997 


0.288 


-124.270 


0.283 


Aur 


54,491.6835 


0.047 


85.742 


0.159 


-123.753 


0.159 


/3 Aur 


54,491.7588 


0.066 


81.639 


0.350 


-119.413 


0.341 


Aur 


54,491.7881 


0.074 


79.167 


0.413 


-116.398 


0.395 


Mizar 


54,488.9333 


0.660 


21.818 


0.361 


-37.796 


0.408 


Mizar 


54,488.9434 


0.661 


22.811 


0.980 


-37.749 


1.093 


Mizar 


54,488.9518 


0.661 


22.457 


0.312 


-37.301 


0.394 


Mizar 


54,488.9604 


0.662 


23.821 


0.556 


-37.302 


0.708 


Mizar 


54,492.0078 


0.810 


44.126 


0.247 


-57.208 


0.329 


Mizar 


54,492.0468 


0.812 


43.908 


0.354 


-57.704 


0.456 


Mizar 


54,548.8658 


0.578 


11.679 


0.482 


-26.290 


0.640 


Mizar 


54,549.8253 


0.625 


17.287 


0.174 


-32.974 


0.222 


Mizar 


54,576.7827 


0.937 


39.127 


0.417 


-52.615 


0.543 


Mizar 


54,578.7569 


0.034 


-74.826 


0.211 


59.789 


0.266 


Mizar 


54,579.7472 


0.082 


-81.159 


0.264 


66.865 


0.333 


Mizar 


54,579.7558 


0.082 


-81.404 


0.289 


66.645 


0.366 


Mizar 


54,602.7619 


0.202 


-52.427 


0.283 


37.381 


0.368 


Mizar 


54,603.6351 


0.245 


-42.830 


0.240 


27.669 


0.301 


Mizar 


54,604.6286 


0.293 


-33.460 


0.289 


19.960 


0.377 


Mizar 


54,634.6759 


0.756 


36.381 


0.306 


-49.740 


0.402 


Mizar 


54,636.6348 


0.851 


48.098 


0.184 


-61.666 


0.246 


Mizar 


54,637.6408 


0.900 


48.821 


0.213 


-62.867 


0.287 


K Peg 


54,400.6227 


0.069 


26.262 


0.031 


2.603 


0.361 


k Peg 


54,401.6788 


0.245 


-10.951 


0.035 


-0.613 


0.345 


k Peg 


54,402.6158 


0.402 


-46.681 


0.055 


4.407 


0.415 


k Peg 


54,402.6305 


0.405 


-47.035 


0.064 


3.843 


0.465 


k Peg 


54,403.6430 


0.574 


-49.989 


0.039 


0.649 


0.309 


k Peg 


54,404.5841 


0.732 


-16.839 


0.042 


-0.683 


0.416 
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Because our template spectra are generated using atomic transition wavelengths from 
the NIST catalog, the derived RVs can be considered "accurate" in the sense that they 
reflect the total Doppler shift between the rest wavelength of a line and the observed 
wavelength. We have not (yet) attempted to tie our velocity scale to any IAU velocity 
standards, nor do we make any correction for gravitational redshift effects. The only 
adjustment made to the RV data is conversion to a solar system barycenter reference frame, 
using the IRAF tool bcvcorr. These barycentric RV data are listed in Table |2j It should 
be noted that for the k Peg RVs, "Vi" refers to the Bb component and "V2" refers to the A 
component. 

We derived orbital parameters from our RV data points using the IDL routines 
CURVEFIT, a gradient-expans ion nonlinear lea st-squares fitting algorithm included with the 



IDL package, and HELI0_RV ((Landsman 



19931 ). which computes a line-of-sight velocity curve 



for a binary component given the period P, periastron time T (or for circular orbits, the time 
of maximum positive velocity), eccentricity e, periastron longitude u, RV semi-amplitude 
K, and systemic velocity Vq (alternatively denoted as 7 by some researchers). We fit the 
primary and secondary RV points simultaneously, assuming that P, T, e, and Vq are the 
same for both components, and that u% and U2 differ by 180°. The CURVEFIT routine 
returns la uncertainties (standard deviations) for all derived parameters. For all of our 
SB2 targets, we adopted the orbital period P from previously-published analyses, because 
our observations covered a relatively short period of time. We did not correct for the light 
travel time across each binary system, because the resulting changes in the RV values are 
small compared to the RV error bars in all six cases. 
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Table 2 — Continued 



star 


HJD - 2,400,000 


phase 


Vi (km s- 1 ) 


V\ error (km s 1 ) 


V 2 (km s- 1 ) 


V2 error (km s 1 ) 


k Peg 


54,404.5988 


0.734 


-16.173 


0.054 


0.520 


0.506 


k Peg 


54,634.8531 


0.294 


-24.846 


0.029 


0.465 


0.317 


k Peg 


54,635.8924 


0.468 


-55.088 


0.032 


-0.740 


0.284 


k Peg 


54,636.8741 


0.632 


-41.896 


0.035 


6.552 


0.275 


k Peg 


54,637.8349 


0.793 


-2.007 


0.034 


0.883 


0.333 


rj Vir 


54,491.9611 


0.229 


11.665 


0.022 


-11.833 


0.067 


rj Vir 


54,548.8113 


0.021 


-27.988 


0.034 


40.160 


0.101 


?7 Vir 


54,549.7519 


0.034 


-26.057 


0.030 


37.392 


0.090 


?7 Vir 


54,576.7011 


0.409 


21.091 


0.035 


-25.264 


0.099 


rj Vir 


54,577.7526 


0.424 


21.138 


0.035 


-25.293 


0.105 


?7 Vir 


54,579.7025 


0.451 


20.923 


0.039 


-25.080 


0.119 


rj Vir 


54,603.7043 


0.785 


-6.225 


0.028 


10.650 


0.085 


?7 Vir 


54,604.6471 


0.799 


-8.312 


0.029 


13.633 


0.087 


77 Vir 


54,605.6656 


0.813 


-10.717 


0.030 


16.643 


0.091 


rj Vir 


54,634.6939 


0.217 


9.759 


0.033 


-10.524 


0.101 


?7 Vir 


54,636.6521 


0.244 


12.931 


0.026 


-14.812 


0.076 
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3. Results on double-lined binary systems 



3.1. i Pegasi 

Our RV measurements for the double-lined spectroscopic binary i Pegasi (HR 8430, HD 
210027, HIP 10917 6) are plotted in Figure [p The most recent published RV work on this 



system comes from iFekel &: Tomkinl (Il983l ). whose orbital parameters are listed in Table [3] 
along with the values that we derive from our dFTS2 observations. In addition to adopting 
their value of the system's orbital period, we also followed their lead in assuming a circular 
orbit, because our RV points only covered half of the orbital phase, and CURVEFIT could 
not place meaningful constraints on e or u. Our values for K\ and K2 are compatible with 
those of Fekel & Tomkin, although our solution for Vq differs by a statistically significant 
amount. This discrepancy may indicate the grav itational influence of a unseen and distant 



Tokovinin et al. 



( 120061 ) did not find any 



third stellar component of the system, although 
close tertiary comp anions in 2MASS images of 1 Peg, and the astrometric observations of 



Bo den et al 



( 119991 ) saw no evidence for a compansion either. Alternatively, the difference 
in Vo might merely be a result of different RV zero points between Fekel & Tomkin's 
observations and ours — unfortunately, we did not make any RV observations of 1 Piscium, 
the RV standard star that they used as their reference spectrum. 

Although the RMS scatter of our RV points around the best-fit orbital curves is small 
(56 m/s and 158 m/s for the A and B components, respectively), the scatter is larger 
than would be expected from the error bars on each individual RV measurement, and is 



significant 



Behretal 



above the instrumental RV error floor of ~ 10 m/s that we determined in 



( 120091 ). To account for this discrepancy, we multiply the RV error bars by 3.11 



for the primary and 1.83 for the secondary, such that the mean per-measurement error 
bar matches the RMS deviation ctrv for each stellar component. The orbital fits are then 
recalculated using these scaled error bars, and the resulting orbital parameters are listed in 



13 



the rightmost column of Table [31 This same procedure is applied to all subsequent binary 
systems as well. 

The additional RV variability, if real, may be a result of stellar activity on both the 
primary and secondary, driven by tidal interactions between the two stars. Fekel & Tomkin 
measure vsini = 7 ± 2 km/s for the primary, very close to an estimated synchronous rate 
of 6.5 km/s, and vsini = 9 ± 3 km/ s for t he secondary, which is well above the estimated 
synchronous rate of 4.5 km/s. iGrayl ( 119841 ) finds vsini values of 6.5 ± 0.3 km/s (primary) 
and 5 ± 1 km/s (secondary), suggesting that the system is synchronized. Our preliminary 
analysis of line broadening indicates vsini = 7.6 and 7.2 for the primary and secondary, 
respectively. If the secondary is indeed spinning more rapidly than the synchronous rate, 
then above-average surface activity could result, which would add significant astrophysical 
RV "jitter" to our measurements. A synchronously-rotating component would be less 
susceptible to tidal effects, but activity migh t still be enhanced by the proximity of a 



massive companion. However, 



Konacki et al. 



( 120091 ) measured RVs of i Peg with three 



different spectrographs, and found no jitter greater than 17 m/s (primary) and 85 m/s 
(secondary), suggesting that the jitter observed by dFTS2 was instrumental rather than 
astrophysical. 



Bodenetal 



(119991 ) measured an inclination angle for this system of % = 95.67° ± 0.22 c 



(based on their primary data set). Using this valu e along with our or bital parameters 



Torres et al 



(120 lOl ). we derive stellar 



with the fundamental parameters recommended by 
masses of M 1 = 1.3241 ± 0.0018 M© and M 2 = 0.8251 ± 0.0010 M , which represent 
relative (statistical) errors of 0.14% and 0.12% respectively. With scaled error bars, the 
mass estimates are the same, albeit with larger error bars, for relative uncertainties of 
0.19% and 0.15%. These values agree reasonably well with the calculations of Boden et al., 
who used Fekel & Tomkin's Ki and K 2 values to determine Mi = 1.326 ± 0.016 M and 
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M.2 = 0.819 ± 0.009 Mq. For our mass estimates, the largest component of the error budget 
is due to the uncertainty in i, although the uncertainties in the K values are also significant 
contributors. 



3.2. uj Draconis 



The spectros copic orbit of the u Dr aconis system (HR 6596 , HP 160922, HI 



was measured by 



Mayor &: Mazehl (119871 ) and more recently by 



Fekel et al 



86201' 



(|2009|). Their 



derived orbital parameters are shown in Table HJ along with our values. Our K velocities 
agree closely with those of Fekel et al. We find a small but nonzero eccentricity for the 
orbits, 0.0023 ± 0.00 02. Fekel et al. derived e = 0.0027 ± 0.0008 for the primary, with 



40.1° ± 17.8° ((Fekel 



20091 ). but their e and u values for the secondary did not 



agree with those of the primary, so they adopted a circular orbit for the system. (Our 
U\ = 137.86 ± 13.48, as described below.) Despite the similarity between our e value and 
their nonzero e value, the measurements of uj\ are substantially different, so we cannot 
plausibly claim that an orbital eccentricity has been clearly detected. 

As an additional test of the nonzero eccentricity, we follow the Fekel et al. procedure 
of fitting orbits to the A and B components separately and comparing the derived 
values for u\ and 0J2, which should differ by 180°. Using the formal error bar data, we 
calculated e x = 0.0020 ± 0.0003 with u x = 74.16° ± 12.74°, and e 2 = 0.0027 ± 0.0005 with 
u>2 = 243.14° ± 13.05°. The eccentricity values agree reasonably well, and the u angles differ 
by ~ 169°, which is within la of 180°. However, this u\ value does not agree with the u\ 
value from the combined fit, which is puzzling. This discrepancy may be related to the 
apparent systematic trends in the secondary RV residuals which are evident in Figure [21 
Further observations with better phase coverage will be required to validate or refute our 
measured eccentricity for u Dra. 
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Table 3. Orbital parameters and stellar mass estimates for t Pegasi. 



parameter 


Fekel & Tomkin (1983) 


this work (formal RV errors) 


this work (scaled RV errors) 


P (days) 


10.213033 ±0.000013 


adopted from F&T 


adopted from F&T 


l ^reaucea rijuj 


40OZU. 14zo 


t)4oyy.ozyD ± u.uuuo 


^/Iqog f;oqq -i- n nnr»7 
O4oyy.ozoo ± u.uuu* 


e 


0.0 assumed 


0.0 assumed 


0.0 assumed 


Kx (km/s) 


48.1 ±0.2 


48.380 ± 0.006 


48.380 ±0.018 


K 2 (km/s) 


77.9 ±0.3 


77.637 ±0.027 


77.638 ±0.050 


V (km/s) 


-5.5 ±0.2 


-4.245 ±0.007 


-4.229 ±0.015 


N ohs 


32 


16 


16 


X 2 primary 




142.20 


16.45 


X 2 secondary 




64.52 


16.11 


ctrv primary (km/s) 


0.90 


0.056 


0.059 


ctrv secondary (km/s) 


1.16 


0.158 


0.156 


M 1 (M ) 


1.326 ±0.016 a 


1.3239 ±0.0018 


1.3239 ±0.0025 


M 2 (M ) 


0.819 ±0.009 a 


0.8250 ±0.0010 


0.8250 ±0.0013 



from 



Bodenetal 



(|1999l ). using F&T velocities in conjunction with spatial interferometer observations 



- 16 - 



As with i Peg, our value for Vq differs from the prior work by a statistically significant 
amount, although the magnitude of the difference is not as large. Differences in the RV zero 
point are the most likely explanation. 

No visual orbit or estimate of the inclination angle i has b een determined for this 



binar y, despite efforts to resolve it using speckle interferometry (jlsobe 



1991 



Miura et al. 



19951 ). We are therefore unable to calculate the true masses of the stellar components. We 
find Mi sin 3 i = 0.16054 ± 0.00011 M and M 2 sin 3 % = 0.13030 ± 0.00007 Af , in moderately 
good agreement with Fekel et al. (With the scaled RV error bars, our mass estimates are 
virtually unchanged, with error bars approximately three times larger.) We hope that 
long-baseline interferometers will soon be able to resolve the astrometric orbit of this system 
and determine the inclination angle. 



3.3. 12 Bootis 



The spectroscopic binary 12 B ootis (HR 53 4, HP 123999, HIP 69226) has 



received recent attention from both 
and astrometric data, and 



Boden et al. 



20051 1. who combined spectroscopic 



Tomkin fc Fekell (120061 ). who performed a high-precision 



spectroscopy-only assessment of the orbit. Orbital parameters are shown in Table and 
our RV data are plotted in Figure [3j The derived quantities for e and U\ are in excellent 
agreement among all three studies. Our K\ and K 2 values, on the other hand, are smaller 
than those of Boden et al. and Tomkin & Fekel by several standard deviations, and our 
derived systemic velocity is different as well. The discrepancy in Vq may simply be ascribed 
to a different RV zero point, but the difference in K\/2 deserves further scrutiny. Due to the 
premature conclusion of our observing program, our RV data do not fully cover the region 
of maximum absolute velocities around phase = 0.15, so the velocity amplitudes are not as 
reliably constrained as they might be. When dFTS observations resume, 12 Boo will be one 
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Table 4. Orbital parameters and stellar mass estimates for u Draconis. 



parameter 



Mayor fc Mazeh (1987 



Fekel et al. (2009 



this work (formal RV errors) this work (scaled RV errors) 



P (days) 

T (reduced HJD) 
e 

OJ! (deg) 
K x (km/s) 
K 2 (km/s) 
V (km/s) 
N ohB 

X 2 primary 

X 2 secondary 

°"RV primary (km/s) 

"RV secondary (km/s) 



5.279799 ± 0.000003 
44698.273 ± 0.005 


35.8 ±0.3 
45.2 ± 0.3 
-14.1 ± 0.2 
27 



5.2798088 ± 0.0000083 
53980.1606 ±0.0006 
0.0 assumed 

36.326 ± 0.029 
44.699 ± 0.039 
-13.975 ± 0.018 
82 



0.19 (unit weight) 



adopted from Fekel et al. 
54547.1347 ±0.0753 
0.0023 ± 0.0002 a 
139.01 ± 5.14 
36.293 ± 0.008 
44.717 ±0.014 
-13.497 ±0.006 
16 
92.10 
138.65 
0.061 
0.126 



adopted from Fekel et al. 
54547.1180 ± 0.1974 
0.0023 ± 0.0006 
137.86 ± 13.48 
36.292 ± 0.020 
44.718 ±0.038 
-13.501 ±0.016 
16 
15.42 
19.42 
0.059 
0.126 



Mi sin 3 i (Mq ) 
M 2 sin 3 i (Mq) 



0.163 ±0.003 
0.129 ±0.002 



0.16090 ± 0.00032 
0.13076 ± 0.00024 



0.16054 ±0.00011 
0.13030 ±0.00007 



0.16054 ±0.00030 
0.13029 ±0.00018 



a but see text regarding the validity of this nonzero eccentricity 
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of our highest-priority targets, so that this issue can be addressed. 

Given smaller K amplitudes than prior publications, we derive smaller masses as well. 
Using i = 107.990° ± 0.077° from Boden et al, we determine M x = 1.4013 ± 0.0025 M and 
M.2 = 1.3607 ± 0.0024 M , for a relative statistical uncertainty of 0.18%. These values are 
several a smaller than the masses derived by Boden et al. and Tomkin & Fekel. 



3.4. V1143 Cygni 



The eclipsing SB2 system of VI 14 3 Cygn i (HR 7484, HD 185912, HIP 96620) was 



previously analyzed by 



Andersen et al. 



(119871 ) . Their orbital parameters are compared 
to ours in Table [6J and our RV data are plotted in Figure HI There is broad agreement 
between the two sets of orbital elements, although the error bars that we derived for K\ 
and K2 are relatively large, partly because of the small number of observations, partly 
because the individual RV measurements had larger error bars due to larger rotational 
broadening of the absorption lines and lower signal-to-noise ratio. Of particular interest is 
the com parison of the periastron ang le u\. From precise photometric timing of the system's 



eclipses, iGimenez &: Margrave! (119851 ) detected apsidal motion (precession of the periastron 
point) with a period of 10,750 years, which would imply a change in the value of oj\ of 0.76° 
during the ~ 22.6 years that elapsed between their last observations (October 1985) and 
our first observations (May 2008). The actual measured change in oj\ is +0.39° ± 0.25°. The 
predicted periastron precession is therefore not ruled out, but is not solidly confirmed either. 
With more extensive observations of VI 143 Cyg, we hope to place more useful constraints 
on the magnitudes of the classical gravitational quadrupole and general relativity effects 



which cause the precession. 
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Table 5. 


Orbital paramet 


ers and stellar mass 


estimates for 12 Bootis. 




parameter 


Boden et al. (2005) a 


Tomkin & Fekel (2006) 


this work (formal RV errors) 


this work (scaled RV errors) 


P (days) 


9.6045492 ± 0.0000076 


9.6045529 ± 0.0000048 


adopted from T&F 


adopted from T&F 


T (reduced HJD) 


51237.7729 ± 0.0051 


52400.4292 ± 0.0035 


54542.2431 ± 0.0031 


54542.2424 ± 0.0042 


e 


0.19233 ±0.00086 


0.19268 ±0.00042 


0.1928 ± 0.0003 


0.1928 ±0.0004 


Wl (deg) 


286.67 ±0.19 


286.87 ±0.14 


286.79 ±0.12 


286.78 ± 0.17 


I<! (km/s) 


67.302 ± 0.087 


67.286 ± 0.037 


67.107 ±0.035 


67.113 ± 0.047 


X 2 (km/s) 


69.36 ±0.10 


69.30 ±0.05 


69.110 ±0.037 


69.102 ± 0.054 


V (km/s) 


9.551 ±0.051 


9.578 ±0.022 


10.040 ±0.018 


10.046 ± 0.025 


•Nobs 


19 


24 


13 


13 


X 2 primary 


~ 49.0 




17.87 


10.16 


X 2 secondary 


~ 49.0 




30.79 


14.02 


°"RV primary (km/s) 


0.47 


0.11 


0.082 


0.077 


orv secondary (km/s) 


0.54 




0.115 


0.119 


Mi (M ) 


1.4160 ± 0.0049 


1.416 ±0.003 


1.4013 ± 0.0025 


1.4011 ±0.0031 


M 2 (Af ) 


1.3740 ± 0.0045 


1.375 ±0.002 


1.3607 ± 0.0024 


1.3608 ±0.0028 



a combined fit to RV and astrometric data 
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Andersen et al 



photo metry by 



987) ad opt an inclinatio n angle of 87 .0° ± 1° based upon ec lipse 



Woodl (I197ll ). 



Popper fc Etzel (I198lh . and 



van Hamme &: Wilson 



(11984] ) . Using that same value for i, we estimate Mi = 1.3815 ± 0.0114 M and 

M.2 = 1.3451 ±0.0100 M Q . The error bars on K\ and Ki dominate the error budget for the 

masses, so further high-accuracy spectroscopic observations of this system are clearly called 



for. 



3.5. (3 Aurigae 



(3 Aurigae (HR 2088, HD 40183, HIP 28360) is an other e c 



consisting of two A2 subgiants in a close 



curves of both components, and 



-day orbit. 



ipsin g double-lined binary, 



Smithl (119481 ) measured the RV 



Pourbabd (120001) reanalyze d these data in conjunction with 



interferometric astrometry data from lHummel et al.l (119951 ) to refine the orbital parameters. 
Our RV curve is displayed in Figure The phase coverage was insufficient to constrain 
e or Wi, so we assumed a circular orbit. Our derived parameters differ from the prior two 
analyses (Table [7]), with K\ and K2 semi-amplitude values intermediate between those of 
Smith and those of Pourbaix. As with VI 143 Cyg, this system exhibits rotational line 
broadening of 30-40 km/s, which increases the uncertainty of each RV measurements and 
thus the derived orbital parameters. 



Adopting i = 76.0° ± 0.4° from 



Hummel et al. 



(119951 ). we calculate M x 



2.3885 ± 0.0134 M and M 2 = 2.3270 ± 0.0130 M Q . (Pourbaix uses i = 75.0° ± 0.73°; the 
source of this value is unclear.) The majority of the 0.54% relative error in our mass values 
is due to the inclination angle uncertainty, so this system would be a prime follow-up target 
for further long-baseline spatial interferometry. 
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Table 6. Orbital parameters and stellar mass estimates for VI 143 Cygni. 



parameter 



Andersen et al. ('1987') this work (formal RV errors) this work (scaled RV errors) 



P (days) 

T (reduced HJD) 

e 

u)x (deg) 
K l (km/s) 
K 2 (km/s) 
V (km/s) 
N ohs 

X 2 primary 

X 2 secondary 

ctrv primary (km/s) 

ctrv secondary (km/s) 



7.64075217 ±0.00000051 
42212.76652 ±0.00015 
0.540 ±0.003 
48.6 ±0.02 
88.20 ±0.20 
91.10 ±0.40 
-16.5 ±0.7 
62 



1.1 

2.2 



adopted from Andersen 
54598.1853 ± 0.0027 
0.5469 ±0.0010 
48.99 ±0.25 
88.867 ±0.248 
91.267 ±0.311 
-16.505 ±0.074 
7 

63.01 
46.49 
0.734 
1.012 



adopted from Andersen 
54598.1835 ±0.0088 
0.5484 ±0.0032 
48.84 ±0.83 
89.055 ± 0.836 
91.508 ±0.941 
-16.461 ±0.235 
7 
5.98 
5.04 
0.789 
0.951 



Mi (AT©) 
M 2 (M ) 



1.391 ±0.016 
1.347 ±0.013 



1.3815 ±0.0114 
1.3451 ±0.0100 



1.3868 ±0.0334 
1.3496 ±0.0308 
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Table 7. Orbital parameters and stellar mass estimates for j3 Aurigae. 



parameter 



Smith (1948) 



Pourbaix (2000 



this work (formal RV errors) this work (scaled RV errors) 



P (days) 

T (reduced HJD) 
e 

Wl (deg) 
K x (km/s) 
K 2 (km/s) 
Vo (km/s) 

A^obs 

X 2 primary 

X 2 secondary 

°TtV primary (km/s) 

(Trv secondary (km/s) 



3.9600421 ± 0.0000013 
31076.719 
0.0 
0.0 

107.46 ±0.39 
111.49 ±0.37 
-17.06 ±0.27 
21 



3.96004 ± 0.00000267 
43915.7 
2.75266 x 10" 6 ± 0.007 
139.043 ± 360.0 
110.246 ± 1 
110.52 ±2.1 
-15.7536 ±0.62 
21 



2.740 
6.369 



adopted from Pourbaix 
54539.0162 ± 0.0003 
0.0 assumed 

108.053 ± 0.072 
110.911 ±0.071 
-17.552 ±0.037 
20 

54.19 

40.83 

0.363 

0.358 



adopted from Pourbaix 
54537.0362 ± 0.0004 
0.0 assumed 

108.053 ± 0.099 
110.911 ±0.098 
-17.552 ±0.052 
20 

28.30 

21.15 

0.363 

0.358 



Mi (M ) 
M 2 ( Mq ) 



2.4 ±0.1 
2.44 ± 0.073 



2.3885 ±0.0129 
2.3270 ±0.0126 



2.3885 ± 0.0134 
2.3270 ±0.0130 



a These orbital parameters are listed in the downloadable data table at the SB9 website (http:/ /sb9. astro. ulb.ac.be/mainform.cgi ) but cannot 
be accessed directly from the web interface. 
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3.6. Mizar A 



The brighter component of the visual binary Mizar (HR 5054, HD 116656, HIP 65378) 



is itself a spectroscopic binary, with two e arly- A dwarfs in an ell i 
lists the orbital param eters measured b y 
revisions computed by 



ptica' 



xj) wi 



Fehrenbach fc Prevot 



Pourbabd (120001 ) with astrometry data from 



20-day orbit. Tabled] 



(119611) and the sub s equen t 



Hummel et al. 



(1998) 



along with the values computed from our RV data (Figure [6]). We find smaller velocity 
amplitudes and a slightly larger eccentricity for this binary system than prior researchers. 
Like the prior two targets, Mizar A's component spectra are moderately rotationally 
broadened, reducing the quality of the RV measurements. 

According to Hummel et al., i = 60.5° ± 0.3° for this binary. Using this value, we 
determine that Mi = 2.2224 ± 0.0221 M & and M 2 = 2.2381 ± 0.0219 M , for a relative 
mass error of 1.00% and 0.98% respectively. The uncertainty in i is responsible for most of 
the mass error; if <jj could be reduced to 0.05° and the errors on K\ and K 2 can be cut in 
half, then the mass uncertainty would drop below 0.25%. 



4. Results on triple systems 



4.1. k Pegasi 



The k Pegasi system (HR 8315, HD 206901, HIP 107354) is a hierarchical triple, with 
two bright components (A & B) in a 11.5-year orbit, and a fainter unseen component in a 
6-da y orbit around th e B co mponent. The canonical published orbit for t his sy stem comes 

(with orbit 



from 



Mavor fc Mazehl ( 



analyses published in 



1987) , with more recent obser vations by 



Mutcrs paugh et al. 



dFTSl prototype (IHajian et al. 



(120061 ) and 



Konackil (12005) 



Muterspaugh et al. 



( 120081 )) and our 



20071 ). Table [9] displays the orbital parameters as measured 



by Mayor Sz Mazeh, dFTSl, and dFTS2. Following Mayor Sz Mazeh, we assume that e = 
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Table 8. Orbital parameters and stellar mass estimates for Mizar A. 



parameter 



Fehrenbach fc Prevot (1961 



Pourbaix (2000 



this work (formal RV errors) this work (scaled RV errors) 



P (days) 

T (reduced HJD) 
e 

OJ! (dcg) 
Ki (km/s) 
K 2 (km/s) 
V (km/s) 
N obB 

X 2 primary 

X 2 secondary 

CT RV primary (km/s) 

"RV secondary (km/s) 



20.5386 
36997.212 ± 0.022 
0.537 ±0.004 
104.16 ± 1.15 
68.80 ± 0.79 
67.60 ±0.91 
-5.64 ±0.15 
15 



1.87 
1.32 



20.5385 ± 0.00013514 
38085.7 ± 0.0269224 
0.529404 ± 0.0052 
105.5 ±0.79 
67.2586 ± 0.96 
69.1796 ±0.77 
-6.3077 ±0.38 
15 



1.88094 
2.39922 



adopted from Pourbaix 
54536.9882 ± 0.0068 
0.5415 ± 0.0010 
105.21 ± 0.14 
66.479 ± 0.095 
66.012 ± 0.118 
-7.342 ± 0.052 
18 
82.55 
81.44 
0.556 
0.641 



adopted from Pourbaix 
54536.9904 ± 0.0106 
0.5415 ±0.0016 
105.27 ±0.23 
66.478 ±0.153 
66.019 ±0.177 
-7.309 ± 0.081 
18 
33.70 
34.01 
0.566 
0.638 



Ah (Mq) 

M 2 (Mq) 



2.5 ± 0.11 
2.5 ± 0.12 



2.2224 ±0.0221 
2.2381 ± 0.0219 



2.2228 ±0.0250 
2.2383 ± 0.0246 



with some values from the SB9 catalog I Pourbaix et al.1 2004) 
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for the short-period orbit. The dFTS2 observations of this system were made during two 
observing runs separated by ~ 0.64 years, so we expect the Vo value of the short-period 
binary to change due to the long-period orbit. To account for this change, we treated the 
data from the two observing runs completely separately, and the results are given in two 
separate columns in the table. For the RV plot in Figure 0, we shifted Bb component's RV 
values from the first observing run by —1.198 km/s so that the two different epochs would 
share the same Vo- 

With only four RV measurements in the second observing run, the value for the Ki 
amplitude should be considered provisional, but the change in Vo is quite clear. We see an 
even larger change in Vo as compared to the mid-1981 observations of Mayor & Mazeh, 
although differences in RV zero point must be considered. Unfortunately, our prior dFTSl 
observations did not yield a value of Vo, because the template spectra were not referenced 
to an absolute wavelength standard. 

RV measurements of the A component are not as precise as those for the B component, 
because its lines are much broader: we estimate vsini = 47.9 km/s for A and 7.1 km/s for 
B. Even taking this fact into account, however, we find a much larger scatter of our RV 
measurements for A than expected from x 2 statistics. One possible expl anation for this 



discrepancy is that the A component is also a close binary, as proposed by 



Beardsley &: King 



( 119761 ). We phased our RV data to their claimed 4.77- day period, but did not fin d any 



( 120061 ) see no 



coherent cyclic pattern in the A component velocities. ( IMuterspaugh et al. 
evidence for a fourth component either.) The binarity of A is a possibility that we might 
explore with future data, but for the time being, this hypothesis is not supported. 
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by 



4.2. r] Virginis 

The n Virginis triple system (HR 4689, HD 107259, HIP 60129) was studied extensively 



Hartkopf et al 



(Il992l ). who combined spectroscopy and speckle interferometry to 



determine the orbits of both th e short-period (72 day) pair and the long-period (13.1 year) 



grouping. 



Hummel et al. 



(120031 ) made additional observations of this system with the NPOI 
interferometer, refining the orbital parameters and determining the inclination angle of the 
close binary orbit. Table [TO] compares their parameters for the short-period Aa-Ab pair 
to our derivation. The parameters are in general agreement, except that the change in Vo 
is larger than previously seen, and is likely due to the gravitational influence of the third 
component of the system. Figure [H] shows the RV curves for the close binary, as measured 
by dFTS2. Note that the residuals for the secondary RV points all lie above the dotted line 
denoting zero residual. This offset may indicate a significant difference in the gravitational 
redshift or convective blueshift between the two component stars, or it may be an effect 
of the third component. Interestingly, Hartkopf et al. see a similar effect, but with the 
opposite sign, in that the Vo that they derive from the Aa RV curve is ~ 0.4 km/s larger 
than Vo from the Ab component. 

Hartkopf et al. made a tentative spectroscopic detection of a blended Mg II 4481 

Table 9. Orbital parameters for the single-lined B component of k Pegasi. 



parameter 


Mavor & Mazeh C1987) 


Haiian et al. (20071 


Muterspaueh et al. C2008) 


this work (Oct 2007) 


this work (June 2008) 


P (days) 


5.97164 ± 0.00006 


adopted from M8zM 


5.9714971 ± 0.0000013 


adopted from M&M 


adopted from M&M 


T (reduced H.ID) 


44801.589 ± 0.015 


53681.86 ± 0.04 


52402.22 ± 0.10 


54400.2125 ± 0.0005 


54633.0989 ± 0.0005 


e 





adopted from M&M 


0.0073 ± 0.0013 


adopted from M&M 


adopted from M&M 


Kx (km/s) 


42.1 ± 0.3 


41.572 ± 0.257 




42.301 ± 0.025 


42.527 ± 0.037 


V (km/s) 


-0.8 ± 0.2 




-9.40 ± 0.22 


-12.135 ± 0.016 


-13.333 ± 0.022 


V A (km/s) 








1.352 ± 0.851 


2.019 ± 1.876 


N ohs 


30 


9 


30 


7 


4 


CT RV B (km/s) 


1.1 


0.990 


0.035 


0.063 


0.093 


itrv A (km/s) 






0.250 


2.084 


3.249 
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feature from the faint tertiary component, which appeared to be rotationally broadened by 
about 160 km/s. Our instrument bandpass does not include this line, so we were unable to 
verify their detection, but we performed a crude three-dimensional cross-correlation using 
two narrow-lined template spectra plus a broad-lined A2V template. We found no evidence 
of a consistent RV solution for the third component. 



According to the observations of 



Hummel et al 



(120031 ) . the inclination angle of rj Vir 



Aa-Ab is 45.5° ± 0.9°. Combining this number with our orbital parameters, we determine 
that Mi = 2.4818 ± 0.1158 M & and M 2 = 1.8745 ± 0.0874 M Q . The inclination uncertainty 
is the dominant contributor to the error budget, although the velocimetry results can 
certainly be improved with more data points. 



5. Summary 

Our results demonstrate that dFTS technology is well-suited to high-accuracy radial 
velocity measurements of double-lined spectroscopic systems. We have determined the 
orbital parameters of six binary systems, matching or improving the published values for 
the masses of the component stars. We also observed two double-lined triple systems, 
providing some constraints on the nature of their stars. 

For our future observational programs for spectroscopic binary stars, we are motivated 
by an assortment of specific scientific goals for which the capabilities of a dFTS are 
particularly applicable: 

1. The most immediate goal is to continue to improve the accuracy of orbital 
parameters of binary systems, particularly the K amplitudes, and thus measure stellar 
masses more accurately. These advancements in spectroscopic capabilities must proceed 
in parallel with better astrometric measurements, as determined by current and future 
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Table 10. Orbital parameters and stellar mass estimates for rj Virginis Aa-Ab. 



parameter 


Hartkopf et al. ( 1992) 


Hummel et al. (2003) 


this work (formal RV errors) 


this work (scaled RV errors) 


P (days) 


71.7919 ±0.0009 


71.7916 ±0.0006 


adopted from Hartkopf et al. 


adopted from Hartkopf et al. 


T (reduced HJD) 


47583.98 ± 0.25 


52321.4 ±0.3 


54403.7295 ± 0.0938 


54403.6116 ± 0.3406 


e 


0.272 ± 0.009 (Aa) 


0.244 ± 0.007 


0.2518 ± 0.0011 


0.2519 ± 0.0040 




0.258 ±0.012 (Ab) 








wi (dog) 


200.9 ± 1.5 


196.9 ±1.8 


197.96 ± 0.48 


197.21 ± 1.74 


Ki (km/s) 


26.67 ±0.20 


same as Hartkopf et al. 


26.532 ±0.054 


26.606 ±0.198 


K 2 (km/s) 


35.58 ±0.31 


same as Hartkopf et al. 


35.128 ±0.081 


35.236 ±0.273 


V (km/s) 


5.24 ±0.19 (Aa) 


4.9 ±0.2 


1.055 ±0.009 


1.118 ±0.033 




4.85 ± 0.32 (Ab) 








N ohs 


50 


same as Hartkopf et al. 


11 


11 


X 2 primary 






292.62 


25.75 


X 2 secondary 






263.24 


28.37 


°RV primary (km/s) 


1.96 


same as Hartkopf et al. 


0.129 


0.150 


orv secondary (km/s) 


4.12 


same as Hartkopf et al. 


0.231 


0.218 


Mi (M ) 


2.34 ±0.2 


2.68 ±0.15 


2.4818 ± 0.1158 


2.5039 ±0.1246 


M 2 (M ) 


1.95 ±0.2 a 


2.04 ±0.10 


1.8745 ± 0.0874 


1.8907 ± 0.0932 



a assumed 
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long-baseline spatial interferometers. 



2. With high-accuracy RV measurements spanning longer periods of time, we will 
be able to detect and quantify secular changes in bina r ies' sy stemic velocities (Vq) due to 



tertiary companions. As discussed by iTokovinin et al.l ( 120061 ). the presence of a tertiary 



companion has significant implications for the formation of close binaries. 

3. Observations of near-circular binary orbits will confirm or refute small nonzero 
eccentricities, thus providing observational validation for theories of tidal circularization 
and the influence of external gravitational perturbations such as Kozai resonances. 

4. For highly elliptical systems like VI 143 Cyg, long-term observing programs can 
measure changes in periastron angle to test theories of apsidal precession due to classical 
and relativistic effects. 

5. Because the instrumental profile of a dFTS is easy to calculate a priori from the 
delay sampling function, we can measure spectral line broadening very accurately, e.g. to 
determine the projected rotational velocities of stellar components and thus shed light on 
tidal spin-up/spin-down mechanisms. 
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Fig. 1. — Radial velocity measurements of the SB2 system i Pegasi. Filled squares show the 
measured RV of the primary component, and open squares indicate the secondary compo- 
nent. Our observing campaign ended before we were able to complete the phase coverage of 
this system. 
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Fig. 2. — Radial velocity measurements of the SB2 system u Draconis. 
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Fig. 3. — Radial velocity measurements of the SB2 system 12 Bootis. 
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Fig. 5. — Radial velocity measurements of the SB2 system (5 Aurigae. 
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Fig. 6. — Radial velocity measurements of the SB2 system Mizar. 
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Fig. 7. — Radial velocity measurements of the double-lined triple system k Pegasi. Filled 
symbols denote the primary of the short-period pair ("Ba"), while open symbols denote the 
"A" component of the long-period orbit. 
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Fig. 8. — Radial velocity measurements of the double-lined triple system rj Virginis. 
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